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ABSTRACT 



Context. BeppoSAX observed GRS 1915+105 on October 2000 with a long pointing lasting about ten days. During this obser- 
vation, the source was mainly in the p class characterized by bursts with a recurrence time of between 40 and 100 s. 
' Aims. We identify five segments in the burst structure and accumulate the average spectra of these segments during each satellite 
^ ' orbit. We present a detailed spectral analysis aimed at determining variations that occur during the burst and understanding 
■ the physical process that produces them. 
Methods. We compare MECS, HPGSPC, and PDS spectra with several models. Under the assumption that a single model is 
able to fit all spectra, we find that the combination of a multi-temperature black-body disk and a hybrid corona is able to give 
a consistent physical explanation of the source behaviour. 
Results. Our measured variations in KT^i, r, KTin, and iii„ appear to be either correlated or anti-correlated with the count 
, rate in the energy range 1.6-10 keV. The strongest variations are detected along the burst segments: almost all parameters 
' ■ exhibit significant variations in the segments that have the highest fluxes (pulse) with the exception of Rin, which varies con- 
" . I ' tinuously and reaches a maximum just before the peak. The flux of the multi-temperature disk strongly increases in the pulse 
^ , and simultaneously the corona contribution is significantly reduced. 

Conclusions. The disk luminosity increases in the pulse and the Rin — Ti„ correlation can be most successfully interpreted in 
term of the slim disk model. In addition, the reduction in the corona luminosity at the bursts might represent the condensation 
of the corona onto the disk. 
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1. Introduction 

The microquasar GRS 1915-f 105 ()Fender fc Bellonill2004l) 
was discovered in 1992 in the hard X-ray band 
(jCastro-Tirado et al.l Il992l ). and later identified with a 
binary systern of a n orbital period of about 33 days 



(jGreiner et al.l 120011) containi ng a black hole whose est i- 
mated mass is 14.0±4.4 ("Harlafti s fc Greiiieil l2004h 
Its distance, evaluated bv lMirabel fc Rodriguez] (jl994l ) 
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and later revised by Fender et al. ( 1999f ). is about 12.5 
kpc, and a source inclination of 66°- 70°is inferred from 
the apparent super luminal motion observed in the radio 



jet (|Mirabel fc Rodriguez|[l993 iFender et al.lll999t ) 



The so urce has been extensively observed in the X- 
rav band (Fender fc Bellonill2004 iRodriguez et~al] 12008 



Ueda et al.l 120101: IVierdavanti et al.l l2010f ). The X-ray 
emission is characterized by strong v ariability with alter - 
nating bright bursts and quiet phases. Belloni et al. ( 2000f ) 
classified the X-ray temporal and spec tral behaviour of 
GRS 1915+105 into 12 classes, and iKlein-Wolt et al 
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( 2002[) and iHannikainen et al. ( 2005 ) identified additional 
sets of variability patterns. As in the case of many black 
hole candidates (BHC), the spectrum can generally be 
fitted with at least two components: a first one, domi- 
nating at soft X-ray energies (<10 keV), interpreted in 
terms of thermal emission, and a second one extend- 
ing up to several hundreds keV that can be modeled 
with a power law. The thermal component, which has 
a higher temperature than in other BHCs, has been in- 
terpreted as emission from an accretion disk surround- 
ing the black hole. According to the most accreditated 
hypothesis, the disk is a thin standard one extending 
up to a few tens of ki lometers from the central source 
(IFender fc Belloni 2004 ) . How ever, the presence of a slim 
disk ( Abramowicz et al. 19881) has also been advocated to 
describe some of t he spectra observed froni this source 
(lUeda et al.l l2009t IVierdavanti et all l2010t lYadav et all 
1993) ■ The phenomenological power-law model is used 
to describe a Comptonization component from an opti- 
cally thick corona above the disk with an electron pop- 
ulation con^ainmg_eithei^ or non-thermal par- 
ticles (jZdziarski et al.l 120011 ) and where energy can be 



also gained from the gas bulk motion (jTitarchuk et al, 
19971) . In addition to the main spectral components. 



lower intensity features have occasi onally been detected: 



a Compton reflec tion from the disk (jZdziarski et al.ll2001 



Ueda et al.ll2009l). an iron hne present only in a few obser 
vations (iFeroci et al.l[l99i : iMartocchia et al.ll2002l [2006; 
Ueda et al.ll2010l) . 

During the bursts, the main spectral parameters of 
the thermal component exhibit significant variations, 
which have been interpreted in term of the emptying 
and refilling of the inner portion of the accretion disk 
(jBelloni et al.lll997l ). In agreement with the phenomenol- 
ogy observed in other black holes, the spectral evolution 
of GRS 1915-1-105 in all classes can been described as 
tra nsitions betw e en th ree separate states, as suggested 
by iBelloni et aP (|2000l ) from the analysis of RXTE ob- 
servations: a state A with a small time variability where 
a strong black-body like component with a temperature 
>1 keV dominates the spectrum, a state B with a disk 
temperature that is higher that in state A and a red-noise 
time variability on a timescale >1 s, and state C where 
the spectrum is dominated by a power law with a spectral 
index between 1.8 and 2.6 and white and red noise time 
variability on a timescale >1 s. 

The BeppoSAX satellite observed GRS 1915-^105 for 
about ten days in October 2000. For almost the entire 
duration of the pointing, the s ource was in the p c lass ac- 
cording to the classification of Belloni et al. ( 2000h . char- 
acterized by quasi-regular bursts with a slow rise and a 
sharp decrease. We described the time e volution of the 



source 'heartbeat' activity in this period in lMassaro et al. 



(20l3, lereafter Paper I), where, on the basis of Fourier 
and wavelet analysis, we defined the regular and irregu- 
lar modes of this class. The former mode is character- 
ized by periodograms with one or two prominent peaks 
and a relatively small variance in the dominant timescale 



in the wavelet scalograms. In the irregular mode, peri- 
odograms exhibit several peaks and wavelet spectra have 
highly fluctuating timescales. Moreover, pulses of regular 
series have a number of peaks (multiplicity) that is gener- 
ally no greater than two, whereas in the bursts of irregular 
series it increases even to more than four. The two modes 
were rather clearly separated in time: the irregular mode 
belongs mainly to the central part of observation (170-375 
ks from the beginning), whereas the regular mode was ob- 
served in the first and the last portions. 

In this paper, we present a time-resolved spectral anal- 
ysis of a large fraction of the data series of Paper I with 
the principal aim of obtaining more information on the 
spectral changes associated with the bursts. We included, 
in addition, results from a short simultaneous RXTE ob- 
ser vation. A sp e ctral analy sis of the p mode was prese nted 
bv lTaam et aD ()l997l ) and lVilhu fc NevalainenI (|l998l ) us- 
ing data from RXTE. Nevertheless, a complete study of 
the burst processes, which is typical of the p class, is not 
yet available (whil e this p aper was under the referee pro- 
cess iNeilsen et al. ( 2011 ) published their results on the 
analysis of p mode spectra observed with RXTE). 

In our analysis, we applied to the data set several spec- 
tral models that have been found to closely fit the energy 
spectral distributions of GRS 1915-f 105 in various vari- 
ability classes; we identified those able to provide a satis- 
factory general description of the spectra, and then stud- 
ied the variations in their main properties. Moreover, we 
present the results of a much shorter simultaneous RXTE 
observation. 

Data and selection procedures are illustrated in Sects. 
2 and 3, the spectral analysis is presented in Sect. 4, and 
the time evolution of the spectral parameters is given in 
Sect. 5. The analysis of RXTE data is described in Sect. 
6 and the physical discussion of the results is provided in 
Sect. 7. 

2. Observation and data reduction 

BeppoSAX observed the GRS 1915-1-105 in 2000 from 
October 20 (MJD 51837.894) to October 29. As in 
Paper I, we used the data obtained with the Medium 
Energy Concentrator Spectrometer (MEGS), operatin g 



in the 1.3-10 keV energy band (jBoella et al.l 119971 ) 



and the Phoswich Detector Syst em (PDS) operatin g 



in the 15-300 keV energy band (jFrontera et al.l Il997l ). 
Furthermore, we considered the data obtained with the 
High Pres sure Gas Scintilla tion Proportional Gounter 
rHPGSPC: lManzo et al.lll997l ). which worked in the nom- 
inal energy band 4-120 keV, and covers the gap between 
the two other instruments. 

Data retrieved from BeppoSAX archive at the ASI 
Science Data Center are organized into observations con- 
taining several satellite orbits with a period of ^96 m. In 
Paper I, we divided the observations into runs, each cor- 
responding to a continuous observing period, and named 
the runs with capital letters and sequential numbers. In 
this paper, we present the spectral analysis performed over 
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runs having exposure times greater than 2000 s in order 
to have, in all instruments, a number of counts sufficient 
to confine the spectral parameters. Table [1] lists these 52 
runs, together with the code of the observation from where 
they were extracted, the time in seconds from the starting 
time of the first observation (20 October 2000 at UT = 
21'' 26™ 55*), the exposure, the rate of each instrument 
in the energy range considered for the spectral analysis 
(1.6-10 keV for the MEGS, 8-30 keV for the HPGSPC 
and 15-150 keV for the PDS), and the number of bursts 
detected in each series. 

The extraction and reduction procedures for the 
MEGS and PDS data were already described in Paper 
I; for HPGSPG, we used the latest distributed version of 
SAXDAS (v2.3.1) and applied the standard selection crite- 
ria. The tool used to extract HPGSPG spectra (hpprod- 
UCTS vS.l.l) produces some failure in the gain correc- 
tion of data from observations 212260011 and 212260016, 
which were consequently excluded from the analysis. 

The MEGS background, estimated from high Galactic 
latitude 'blank' fields accumulated in the same region of 
the detectors, is neglected, as it has a count rate of a 
level 2.6x10^^ cs^^, which is much lower than the typical 
source flux, even in the faintest states. Both the PDS and 
HPGSPG background were obtained from simultaneous 
measurements thanks to the standard operational rocking 
mode. 

Our spectral analysis was performed with XSPEG v. 11 
using the response matrices, that are available from the 
5eppoSAX archival- We restricted the spectral analysis of 
HPGSPG data to the energy range 8-30 keV where the 
response matrix is more accurate. The MEGS-PDS inter- 
calibration factor was fixed to 0.85, which is its average 
value in the range 0.83-0.89 for point sources of high count 
rate, and the MEGS-HPGSPC factor was limited to vary 
within the range 0.9-1.1 to take into account the discrep- 
ancy in the flux measured by the two instruments. We 
found that HPGSPG measured a flux higher than 

the MEGS. 

The errors in the values of the best-flt parameters re- 
ported throughout the paper and in the tables correspond 
to a 1 cr confldence level for a single parameter, unless 
explicitly stated. 

3. Burst "anatomy" and count rate variations 

In Paper I, we distinguished three main phase intervals 
or segments of the burst profile and investigated the time 
evolution of the count rates for each of them. The most 
prominent feature is the pulse, which generally consist 
of several short peaks. Prior to the pulse, there is the 
slow leading trail (hereafter SLT), which is also called the 
'shoulder' in the literature, during which the count rate 
increases much slower than in the pulse. There is no clear 
separation between these two components, and we took 
it at the time where the count rate reached the level of 

^ http:/ /www.asdc. asi.it/bepposax 



the FWHM of the pulse, which was evaluated by means 
of a Gaussian best-fit. Between the end of the pulse and 
the beginning of the SLT, the count rate reaches its low- 
est level that can be considered as a baseline upon which 
bursts are superposed. We defined a time interval of three 
seconds where it is possible to estimate this baseline level 
(hereafter BL corresponding to the first three seconds of 
the red segment of the light curve shown in Fig. [Ij and 
found that BL remained remarkably stable in each indi- 
vidual time series (see Paper I for more details). In the 
spectral analysis, we essentially followed this division and 
introduced a finer partition. 

The analysis of the BL spectra gave parameters in 
agreement, within the errors, with those of the SLT ones, 
and owing to the lower quality statistics, we did not con- 
sider BL as a separate segment. The signal in the entire 
SLT is strong enough to divide this component into two 
parts of equal length, hereafter called SLT-1 and SLT-2, 
respectively, which are plotted using data points of differ- 
ent colors in Fig.[Tl where a short segment of the series A3 
is reported as an example. The same colors are also used 
in the other figures to permit an easier identification of the 
various time segments. We also divided the pulse into two 
equal time intervals, namely P-1 and P-2, which are indi- 
cated in Fig. 1 and in subsequent figures by blue and cyan 
data points, respectively. An even finer segmentation to 
enable us to analyze individual peaks was not possible be- 
cause the peak shape and duration is highly variable from 
burst to burst. Finally, the ending part of the pulse to 
the initial minimum of SLT was taken as a fifth segment, 
marked in Fig. 1 with light grey points, and labeled FDT 
corresponding to Final Decaying Trail. The starting and 
ending times of these parts in each burst were evaluated 
by means of the same procedure described in Paper I. For 
each data series of MEGS, HPGSPG (whenever available), 
and PDS, we accumulated five spectra, each of them cor- 
responding to the integral of a selected part of the burst. 

The variations in the count rates of the various com- 
ponents in the course of this long pointing were analysed 
in detail in Paper I. However, to take into account the 
finer segmentation of the bursts' profile introduced during 
the spectral analysis, we again tracked the behaviour of 
each part to highlight any possible differences from the 
results presented in Paper I. We considered the same time 
segmentation of the whole dataset in three intervals intro- 
duced in Paper I: interval I from the start time to 1.7x 10'°^ 
s, interval II from this time to 3.8x10^ s, and the last in- 
terval HI is from 4.0x10^ s to 6.0x10^ s. In particular, 
the second interval is characterized by the presence of the 
irregular mode detected in the timing analysis (Paper I), 
and the last interval corresponds to the final high state 
after the rapid increase in count rate detected between 
3.8x10^ s and 4.0x10^ s. Figures [2] and [3] show the aver- 
age rate for the five portions of the bursts for each series 
in the 1.6-10 keV (MEGS) and 15-100 keV (PDS) energy 
ranges, respectively. 

As already discussed in Paper I, the evolution of the 
count rates is not the same for the different parts of the 
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Fig. 1. A snap shot of the light curve in the entire MECS 
energy range of the series A3. Red and green data points 
indicate the SLT-1 andSLT-2 intervals, blue and cyan 
points are for P-1 and P-2, and grey points indicate the 
ending part of the burst (FDT). 



bursts. The behaviours of SLT-1, SLT-2, and FDT are 
similar in the MECS and PDS data sets, and their count 
rates undergo a significant increase between 380 ks and 
400 ks with nearly stable levels before and after this tran- 
sition, as was also found for the BL level (see the lower 
panel of Fig. 15 in Paper I). 

In the MECS data, the count rates of P-1 and P-2 
remain quite stable during the entire length of the obser- 
vation, without any indication of a step increase. We note 
that the relative pulse contribution in the MECS range is 
higher in the first 400 ks of observation (before the step 
increase) as shown in Fig. 15 of Paper I, where we present 
the difference of the rate with respect to the BL level. In 
the PDS data, the rate variation among the five segments 
is lower than in the MECS (the largest difference is 1.7 
compared to 2.5 in the MECS), and the presence of the 
peak contribution is only evident in P-2, for which there is 
a strong increase in the rate during the first 400 ks with a 
bump in the central part of the pointing where the source 
was in an irregular mode. The hardening of P-2 during the 
irregular mode was already reported in Paper I, where we 
noted that PDS peaks were clearly apparent only in the 
final part of the pulse, while in the regular series they were 
not so high. In summary, there is a strong correlation be- 
tween the MECS and the PDS rate, which weakens at the 
pulse, where the MECS rate is almost constant. 

The spectral variations within the five segments cor- 
respond to c yclica l changes in the state defined by 



Belloni et al 



(I2OOOI) : the source is mainly in state C in 
the SLT, moves to A in P-1, then to B in P-2, and, fi- 
nally, returns to C in the FDT section. 

4. Spectral analysis 

As simple spectral models had been shown not to give ac- 
ceptable fits over an energy range from a fraction of a keV 
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Fig. 2. Time history of MECS (1.6-10 keV) count rates for 
the five segments identified along the bursts. Top to bot- 
tom: SLT-1 (red), SLT-2 (green), P-1 (blue), P-2 (cyan), 
and FDT (grey). Vertical dashed lines mark the main in- 
tervals of the observation. 

to about one hundred keV, we considered composite mod- 
els. The majority of the models that we studied had indeed 
been previously used to investigate this source, although 
for different brightness states and variability classes. These 
models include the emission from a standard optically 
thick accretion disk to repro duce the spectrum up to about 
15 keV (diskbb in XSPEC iMitsuda et al]|l984[) . and, at 
higher energies, they consider other spectral functions that 
generally have a significant contribution also at lower en- 
ergies. We tried several models (listed in Appendix A) that 
interpret the expected spectral shape either with a combi- 
nation of phenomenological functions or with laws based 
on theoretical descriptions of a corona- accretion disk sys- 
tem. 

We required that the same spectral model was able to 
fit all the 52 spectra of each segment by changing only the 
values of the parameters. To test the consistency of the 
model with the data, we compared the frequency distri- 
bution of the reduced with the expected parent distri- 
bution. This was evaluated from 10000 simulated distri- 
butions each generated by randomizing a set of 52 val- 
ues from the functions corresponding to the degree 
of freedom of each observed spectrum. We found that a 
model closely fitting all spectra has an expected distribu- 
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Fig. 3. Time history of PDS (15-100 keV) count rates for 
the five segments identified along the bursts. Top to bot- 
tom: SLT-1 (red), SLT-2 (green), P-1 (blue), P-2 (cyan), 
and FDT (grey). Vertical dashed lines mark the main in- 
tervals of the observation. 



tion of the reduced given by a Gaussian-like function 
centered on 1.00±0.02 that has a standard deviation of 
(T=0.10±0.01. As a result of our fits, we found that the re- 
duced distribution has the expected width (O.lOiO.Ol), 
but the mean (1.05±0.01) is shifted even if compatible 
with the expected one to within two standard deviations. 
This shift is mainly caused by the highest quality data 
statistics (such as those of the pulse) and may be pro- 
duced by residual systematics originating either from local 
detector features at the maximum level of 5% (such as the 
4.78 keV residual in the MEGS spectra) or from the uncer- 
tainty in the intercalibration factors (~3%) fixed for the 
PDS to 0.85 (see Sect. 20. This effect cannot be properly 
treated within XSPEG, where systematics are introduced 
as a single factor over the entire fitting range. To take 
this effect into account, we compared our distribution 
with a Gaussian centered on 1.06 and the compatibility 
of the resulting distributions with the expected one was 
confirmed by means of a Kolmologov-Smirnov test. 

We found that the "good model" that most closely 
describes the source behaviour during the burst and pro- 
vides a consistent physical interpretation is the standard 
disk plus a corona that contains an hybrid population of 
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Fig. 4. Frequency histograms of the reduced x for all 
spectra of the "good model" (continuous orange line) 
and of model iv described in Appendix A (dashed line) 
compared with the expected distribution evaluated by a 
MonteGarlo simulation (black line). 



electrons including both thermal and non-thermal distri- 
butions ( Poutanen fc Svensson 19961) . The code relative 
to this model implemented in XSPEG (cOMPPS) includes, 
together with emission from the multi-temperature disk 
and corona, the Gompton reflection from gas whose tem- 
perature and ionization status can be set as fitting pa- 
rameters. It provides acceptable best-fits to the spectra 
of all data sets and the reduced distribution is com- 
patible with the expected one within only one standard 
deviation as shown in Fig. |4l where the COMPPS x^ dis- 
tribution is compared with the one relative to the model 
iv (top panel) described in Appendix A. The model v in 
Appendix A (oiSKBB-f pow-|-bbody) also provides fully 
acceptable best-fits, without any prominent features in the 
residuals, for the spectra of all data sets. However, as ex- 
plained in Appendix A, no simple physical interpretation 
can be given for the black-body component, and therefore 
it was not considered in the subsequent analysis. 

In the COMPPS fits, the electron distribution is as- 
sumed to be thermal below a Lorenz factor above 
which energies are distributed with a power law spec- 
trum of index F up to a Lorenz factor 7^02; fixed at 10^ 
([Zdziarski et al.l 120051 ). The values of 7mm and F, which 
were free to vary during first sets of fitting, were found to 
maintain almost constant values (7mm: mean 1.3, rms 0.7; 
F: mean 4.6, rms 1.6) in all the observations and in the 
various segments and were then fixed to their average val- 
ues. The 7m,ira c l osely agrees with the results obtained by 
Zdziarski et all (|2005l ) for the analysis of RXTE-I-OSSE 
observations, while, F is fixed to a value slightly larger or 
margin ally compatible with t hose from previou s measure- 
ments (jZdziarski et al.ll2005HUeda et al.ll2010l) . 

The disk, surrounded by a spherical corona (param 
geom={) in XSPEG), extends to between 10 and 
1000 Schwarzschild radii with an inclination of 0=70°. 
Gompton reflection from a non-rotating disk (param 
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BetorlO^-lO in XSPEC) with a gas temperature of 10^ 
K and a low level of ionization (param ^=0) is also taken 
into account. 

A reduction in the free parameters led to a broader 
distribution that was still compatible with the ex- 
pected one (see Fig. 2]). The free- fitting parameters were 
thus: the absorbing column density (-/Vh), the tempera- 
ture (kTei), and the optical depth (r) of the corona, then 
the temperature {kTin) and the apparent inner radius 
(Rin) of the accretion disk, and the solid angle of the 
cold reflecting matter to the illuminating source (il) ex- 
pressed in units of 2tt. The values of Rin are computed 
from the normalization of the model, K, using the rela- 
tion K = (i?i„[km]/I?[10kpc])^ x cos{9), where D is the 
distance (12.5 kpc). They are practically equivalent (to 
within a few p ercent) to the ones con i puted using the for- 
mula given by iKubota fc Makishima that, assum- 
ing isotropic radiation, takes into account the conservation 
of the photon number in the Comptonization process. 

The spectral fittings and residuals for the five segments 
of the series E7 are shown in the left panels of Fig. [S] their 
relative models plus components are plotted in the right 
panels of the same figure. 

4.1. The iron line 

The detection of a fluorescence Ka iron line from 
GRS 19 15-1-105 has bee n reported on sever a l oc- 

12002; 



sorber containing materials in grains whose char acteristics 



casions (Martocchia et al 



Martocchia et all l2006t 



, Z dziarski et al, 
Ueda et al.i l2llToli 



2005 



cently, by iTitarchuk fc Seifinal (l2009l) 



and, more re- 

., , , ., . .- In particular, 

Martocchia et al. ( 20021) from the detection of a strong 
line found evidence of rclativistic effects in a black hole 
environment. 

We searched for the presence of the Ka iron line by 
including in the adopted model a Gaussian profile, whose 
central energy was allowed to vary across the range be- 
tween 6.3 keV and 7.1 keV and whose width was fixed to 
0.1 keV. No significant line was observed in any of the an- 
alyzed spectra with a maximum detection of 2.5 a only in 
a few spectra and a 3a upper limit of ~50 eV to the EW. 

4.2. The low energy absorption 

The low energy ab sorption of GRS 1 9 15-1-105 is com- 



plex and variable (Belloni et al.l 200C 


; Lee et al. 


Martocchia et al. 20061 Yadav 20061: 


Ueda et al. 


2OIOI: van Oers et al. I2OI0I). The values of the 



2002; 



2009 



equiv- 

A^H obtained from the analysis 
between 2x10^^ and 16x10^2 



alent column density 

of X- ray data ranges 

cm-^ dGreiner et al.lll996l:lTaam et al.lll997l:lTV"udolYuboy 



2001: 'M artocchia et al.H2002t llee et alj|2002nYadavll2006t 
van Oers et al.l 2010l) and are generally higher than the 
values derived from the 21 cm hydrogen l ine emission 
along the line of sight f4.7ib0. 2xl0^^ cm"^: IChaty et al, 
19961) (> 1.8 X 10^2 cm-2; iDickev fc LockmanI \l99^ . 



have been derived with Chandra observations (jLee et al 
2OO2I : lUeda et al.ll2009t ). 



These results indicate that there is a composite local ab- 



Since our data are relative to an instrument of mod- 
erate energy resolution, and a low energy boundary at 
~2 keV, we adopted the XSPEC wa bs function assum- 
i ng st andard solar abundance from Anders fc Grevess3 
()l989l) . No significant variation in the absorbing col- 
umn was discerned during the ten days of observations. 
However, we found a rather weak increase in the column 
density at the peak: the A^h varied between 4.3 x 10^^ cm"^ 
in the SLT and 5.4x10^^ cm"^ in the pulse (see Table[2]). 
Fixing the A^h to an intermediate value resulted only in 
a small worsening of the distribution without affecting 
the values of the other parameters. In particular, fixing 
to 5.0x10^2 cm~2 the A^h of the P-2 spectra, differences 
smaller than 10% were measured for the disk parame- 
ters and values that were compatible to within 2a for the 
corona ones. 

Following lUeda et al.l ()2010l ). we evaluated the contri- 
bution of the halo produced by dust scattering to the 
MEGS spectra below 3 keV. We found that it produces 
about 10% of the emission ass uming the fractio nal in- 
tensity derived for GX 13+1 (|Smith et al.ll2002l) . This 
could account for the apparent variability of A^h on a 
long timescale, but the scattered flux is not expected to 
vary during the period 40-100 s coherently with the burst, 
where the time delay respect to the source emission is 
about 100 hours (|Ueda et al.ll2010t) . 

We also fitted our spectra with a more complex absorp- 
tion model (tbvarabs in XSPEC ) by fixing the a bun- 
dances to the va lues obtained by iLee et ah (2002) and 
Ueda et al. I (l2009l) . and in both cases a broader distri- 
bution centered at ~1.2 was obtained. Leaving the abun- 
dances free to vary produced poorly constrained values. 
No clear explanation can be found for this result that may 
also be due to the time variability of the absorber. 



5. Evolution of the spectral parameters 

We analyzed the variations in the best-fit spectral parame- 
ters in the course of the observation to search for possible 
correlations between them and with the changes in the 
mean source luminosity. We present our results in term of 
the two main components of the model: the disk and the 
corona. 

For the three intervals we evaluated the average val- 
ues of the parameters of each component and their rms, as 
given in Table[2j to take into account the variability, which 
is also listed in the table as uncertainties. The time evolu- 
tion of the parameters for the corona and the disk compo- 
nent are shown in Figs.jSland [71 respectively. Moreover, at 
the end of the fitting procedure, we derived the disk con- 
tribution to the total flux by setting r and fl to zero, and 
equated the contribution of the corona to the difference. 
Non-thermal electrons are responsible for ^30% in SLTs 
and 15% in the pulse of the coronal emission as evaluated 
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by setting the parameter gmin to a negative value (i.e., 
purely thermal case). 

The reflection component was detected with a signif- 
icance higher than Str, in the SLT, while no reflection 
is observed in the pulse. Its contribution to the total 
flux is always lower than 10%. The solid angle has al- 
most no signiflcant variations either with time or segment. 
In particular, we obtained 17=1. 2±0. 3 and 17=1. 6±0. 3 
for SLT-1 and SLT-2, respectively, for the average and 
rms value during the ten days of observations. These 
values are higher than expected based on the upper 
limit to the equivalent width of the narrow iron-K line 
(17 <0.5), but, when this parameter is constrained to val- 
ues lower then the upper limit, the maximum (17=0.5) 
for all spectra is found. Several different values are given 
in the literatur e for the reflectin g solid angle, ranging 
from 0.1 to 2 (iFeroci et al.l Il999|: IZdziarski et all l200 



Martocchia et all 120061: lUeda et all 120091: Ivan Oers etal.. 
20inh . that could be related to t he spectral state of the 
source (ISobolewska fc Zvckilliool . Results obtained from 
our analysis are compatib le wi th the ones presented by 
Sobolewska fc Zvckil (|2003h and lZdziarski et al.l (|2005l) for 
other variability classes. However, the low iron abundance 
that we used for our fitting could lead to an overestimate 
of the intensity of the reflection component. 

The variations in kT^i, r, fcTi„, and i?i„ were found 
to be either correlated or anti-correlated with the 1.6-10 
keV MECS and 15-100 keV PDS rates for the majority 
of the segments. These parameters have almost constant 
values in the intervals I and III (regular mode) with a 
significant change in the mean values occurring in cor- 
respondence with the increase in the count rate. During 
interval II (irregular mode), the values are remarkably 
similar to those of interval I, but have a larger variance. 

The strongest variations are detected for the segments, 
as shown in Figs. |5| and HI where the averages for the 
three intervals are plotted with different line styles (a 
solid line for interval I, a dashed line for II, and a dot- 
ted line for III). The figures also show the 0.01-200 keV 
flux in units of 10~^ erg cm~^ s~^ in the bottom pan- 
els computed by extrapolating to lower energy the best-fit 
model. The parameters have very similar values in SLT-1 
and SLT-2, while they exhibit significant variations in the 
pulse with the exception of Rim which varies throughout 
all segments and r eaches a maximum i n SLT -2 as pre- 
viously observed by Vilhu fc NevalainenI ( IQQSl) . The flux 
of the multi-temperature disk strongly increases during 
the pulse, and simultaneously the corona contribution re- 
duces significantly with respect to the values detected in 
SLT. Moreover, the corona increases its emission by 40- 
50% from interval II to interval III, while smaller varia- 
tions (20-30%) are observed in the disk flux during the 
three intervals, suggesting that the increase in the rate 
after 4.0x10^ s is caused mainly by this component. 

The total luminosity of the source, for a distance of 
12.5 kpc, was also computed in the 0.01-200 keV range 
as Ltot = Ldisk + Lcorona, assumiug isotropic emission for 
both components and deriving Ldisk by setting both r and 
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Fig. 6. Corona: Temperature (top figure) and optical 
depth (bottom figure) of the thermal component of the 
electron population in the corona vs time in the five seg- 
ments. 

17 to 0. The source luminosity changes from ~50% of the 
Eddington luminosity (for a black hole mass of 14 Modot) 
to 100% (or slightly above it) in the pulse as shown in the 
second panel of Fig. (TO] In the third and fourth panels of 
the same figure, the fraction of the total luminosity versus 
(vs) segment for the disk and the corona are also shown. It 
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Fig. 7. Multi-temperature black-body disk: Disk 
temperature (top figure) and apparent inner radius (bot- 
tom figure) vs time in tlie five segments. 



is evident that the increase in the pulse flux is due to the 
disk contribution, and that a higher percentage of photons 
is Comptonized by the corona after 400 ks. 
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Fig. 8. Corona: Variability of the spectral parameters 
along the segments. The electron temperature is shown 
in the second panel, the optical depth in the third, and 
the 0.01-100 keV flux in the bottom one. Interval I is rep- 
resented with a solid line, II with a dashed line, and III 
with a dotted line. The top panel shows the MECS rate 
for comparison; two bursts are shown for clarity. 



6. RXTE simultaneous observation 

RXTE observed GRS 1915-f 105 a few times during the 
BeppoSAX long observation. We extracted the data from 
observation 50405-01-04-00, made on 2000 October 27 
06:28-09:25 UT. The starting time of this data series was 
550830.0 s after the beginning of the BeppoSAX. point- 
ing. The exposure times and count rates were 2480 s and 
3033.7 ct s-i for the PGA, and 805 s and 28.3 ct s^^ for 
HEXTE, respectively. Spectral data were available only 
with time resolution of 16 seconds, we then defined from 
the PGA light curve (see bottom panel of Fig. [TT|) two in- 
tervals, SLT (9000-15000 ct s"! ) and P (>15000 ct s'^), 
where were as close as possible to the BeppoSAX finer se- 
lection. Two spectra were extracted from both the PGA 
(PGU2 only) and HEXTE instruments, following standard 
RXTE procedures and a 1.6% systematic error was added 
to the PGA data. Table [5] shows the best-fit parameters 
relative to the "good model" that are in good agreement 
with the BeppoSAX ones. We note that the high-resolution 
light curve, plotted in the to p panel of Fi g. [TJ shows a 



deep "dent" in the bursts fsee lBelloni et al.ll2000l) . This is 



mm . 



true for all RXTE bursts observed during the BeppoSAX 
interval, but since no spectral information was available 
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Phase 

Fig. 9. Multi-temperature black-body disk: 

Variability of the spectral parameters along the seg- 
ments. The temperature is shown in the second panel, the 
apparent inner radius in the third and the the 0.01-100 
keV flux calculated by setting r and J7 to in the bottom 
one. Interval I is represented with a solid line, II with 
a dashed line and III with a dotted line. The top panel 
shows the MECS rate for comparison; two bursts are 
shown for clarity. 
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Fig. 10. Evolution of the total luminosity in units of 
Eddington luminosity (second panel) for the considered 
burst segments. The fraction of the total luminosity for 
the disk and the corona are shown in the third and fourth 
panels, respectively. Interval I is represented with a solid 
line, II with a dashed line, and III with a dotted line. 
The top panel shows the MECS rate for comparison; two 
bursts are shown for clarity. 



for these data, it was impossible to perform a more time- 
resolved analysis. 

7. Discussion 

The phenomenological model used to fit all spectra, has 
allowed us to derive some physical information about the 
characteristics of the source emission in a p mode. The 
scenario that we adopt in the following discussion is that 
of an accretion disk plus a surrounding high temperature 
corona. These two components and their strong interac- 
tion with each other have been proposed by several au- 



thors (iTaam et al 
2004t 



Ueda et all 



1997: Zdziarski et al. 2001: Done et al 



range 1-2 keV and those of Rm are in the range 20- 
70 km. The se results compare quite well with other re- 
sults for o dVilhu fc Nevalainen 1998 ) and other classes 
I Belloni et al.lll997n . althoug h temperatures as low a s 0.7 
keV have also been reported (jFender fc Bellonil 120041 ). 

According to the Shakura-Sunyaev thin disk model, 
Rin depends only on the mass and spin of the black hole. 
Our results are compatible with the gravitational radius 
{Rg — GM/c^) of a 14 Mq black hole, which i s equa l to 
~21 km. As pointed out bv IVierdavanti et"aL I (|2010[) . to 
evaluate the true inner radius we must consider the 
relativistic and s pectral hardening co rrections, which can 
be factorized as (|Kubota et al.lll998l ) 



2009) in describing the complex be- n* _ ^-Sad 



haviour of GRS 1915-^105 



7.1. The multi-temperature disk parameters 

The apparent inner radius and temperature are the two 
most important parameters of the disk component when 
attempting to interpret the emission of GRS 1915-1-105 . 
The values of A;Tj„ found in our analysis vary within the 



(1) 



adopting their values (k = 1.7, C — 0.412), we obtain 
values larger than Rg within the reported uncertainties. 

In any case, our lower values of are smaller than 
the Schwarzschild radius Rs = 2_Rg. This might be an 
indication of either a rapidly rotating black hole or a 
slim accretion disk at about the Eddington luminosity. 
Indications of a spinning black hole in GRS 1915-1-105 have 



10 



T. Mineo et al.: The complex behaviour of GRS 1915+105 in the p class. 



< 
u 




2250 



2300 



2350 
T-Tstart (s) 



2400 



2450 



Fig. 11. Three bursts as observed by the RXTE/PCA on 
2000 October 27 vs time m s from Tstart. The top panel 
shows the light curve with 0.125 s time resolution and the 
bottom panel the 16 s one. Spectral data were only avail- 
able for the 16 s resolution, and the vertical dashed lines 
in the bottom panel indicate the two selection intervals 
{SLT and P) adopted for the spectral accumulations, as 
explained in the text. 



been presented by several authors (jMiddleton et al. I l2006l 
and references therein), who estimated a dimensionless 
spin parameter a = Jc/GM^, where J is the angular 
momentum of the black hole, equal to about ~0.7. With 
this value, a minimum distance of ^36 km is accept- 
able, and considering the large uncertainty in the mass, 
which might be slightly overestimated, our results are 
not in conflict with the general scenario adopted for this 
peculiar sourc e . On the other hand, the calculations of 
Watarai et al. (|2000h for the emission from a slim disk, 
demonstrated that radiation can be produced at distances 
much smaller than the last stable orbit, i.e. in the region 
between Rs and 3Rs- 

The plot of Tin vs for the multi-temperature disk 
obtained for all the fitted spectra is shown in the top panel 
of Fig. [12] The points of SLT-1 and SLT-2 form two clus- 
ters with the same temperature range but clearly sepa- 
rated in terms of Rin, because of the different flux lev- 
els: both clusters exhibit an anticorrelation between the 
two parameters. The tight correlation is clearly evident 
between P-1 and P-2, and appears to behave coherently 
during the entire observation. We fitted the parameter val- 
ues for these two sets with a single power-law 



leaving the exponent g as a free parameter and obtained 
q =0.46±0.01; the fitting relation is shown in Fig. [T2]with 
a solid line. This result translates into an almost constant 
flux at the peak, as shown in the bottom panel of Fig. [T^ 
where the disk luminosity for all the spectra is plotted vs 
hT- 




R. (km) 




CRl 



(2) 



1.5 

KJ,„ (keV) 

Fig. 12. Upper panel: temperature of the multi-color disk 
vs apparent inner radius. Lower panel: disk luminosity, in 
units of Eddington luminosity, vs the inner temperature. 
The continuous black lines are the best-fit models (see 
text). 



Observations of several BHCs have shown that when- 
ever changes have been det ected, the lumino sity of the disk 
varies m proportion to (jPone et al.ll2007i and reference 
within). However, deviations from this relation of the form 
L (X have been observed in some objects, as in the 
case of XTE J1550-564, which has an inclination angle of 
~70° and a mass 8-11 M(d, and exhibits superlumin al jets 
(|Kubota k Dondl2004l : iKubota fc Makishimall2004l ). sim- 
ilar to GRS 1915-1-105. Fitting points from all segments, 
except for that containing the pulse, with a power-law, the 
resulting values for the temperature exponent are 2.7±0.4 
for SLT-1, 2.0±0.5 for SLT-2, and 1.9±0.2 for FDT, re- 
spectively. Several models have been proposed to explain 
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the relation L oa for a slim disk where advection dom- 
inates the energy transfe r in the inner region and reduces 
the radiation efficiency ( Watarai et al. 20001 ) . Moreover, 
the effects of the large inclination angle cannot be ne- 
glected since deviations from the L oc relation are 



main ly observed in highly inclined sources (jPone et al 
l2007h . 



7.2. The corona 

The emission of the corona has been detected up to 
600 keV without any break and can be interpreted 
in terms of both a thermal and a non-thermal elec- 



tron population ([Zdziarski et al.ll2001r) . The output spec- 
trum of the thermal component, which is mainly re- 
sponsible for the spectrum in the BeppoSAX band, de- 
pends on the Comptoniza tion parameter y, defined as 
i|Rvbicki fc Lightmanl[l979h 



y = 4- 



kTpi 



-Max(T, T^), 



(3) 



where me is the electron mass and c the speed of light. 
The values of y vs segment are plotted in Fig. [T31 no large 
variations are observed between the time intervals I and II 
(continuous and dashed lines), while interval III is charac- 
terized by higher values, in agreement with the results for 
the absolute flux emitted by the corona. Significant varia- 
tions in y are observed for the segments with higher values 
in both SLT and FDT, while these variations are smaller 
than 0.1 in the pulse, in accordance with the trend in the 
flux. 

As shown by ISunvaev fc T itarchukI (jl980l) , the ratio 
of the luminosity of the Comptonizing corona Lc to the 
luminosity of the source of the seed photons (Lq) for kT^ > 
hv is a function of the spectral index, thus for two different 
segments we have 



(4) 



^ = ^x/(r„r,), 

where i and j indicate different segments and F is the 
energy index. The spectral index values obtained for the 
model V, are practically constant with F ~1.6 (see model 
V in Appendix A), hence the function f{Tj, F^) is close to 
unity. The decrease in the corona luminosity observed at 
the peak section could then be attributed to a reduction 
in the size of the corona itself. Models of disk+corona sys- 
tems have shown that strong interactions between the two 
components are able to explain the constant spectral in- 
dex in term of a tight interaction between the temperature 
and t he optical depth of the corona ( Haardt fc Maraschi 
19911 ). Moreover, the thermal coupling of the disk and the 
corona would result in either the evaporation of the disk or 
the c ondensation of the corona depending on the accre tion 
rate (|Liu et al.ll201ll: iMever et al.ll2007l:lHaardtl[l993h . 



8. Summary and conclusion 

The spectra of GRS1915-I-105 that we have accumulated 
in five segments of the burst of the p class have been mod- 
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Fig. 13. Corona: Variability of the Comptonization pa- 
rameters y with the burst segments. Interval I is repre- 
sented with a solid line, II with a dashed line, and III 
with a dotted line. The top panel shows the MECS rate 
for comparison; two bursts are shown for clarity. 

eled with several combinations of spectral laws. Only two 
models were found to provide acceptable fits to all spectra: 
a multi-color disk with both a black-body and power law, 
and a disk plus an hybrid corona. Since the former does 
not have a consistent physical interpretation, we studied 
the spectral changes in the source based on the basis of 
the latter. 

The following conclusions can be inferred from the 
analysis of the best-flt values of the parameters: 

— The total 0.01-200 keV luminosity of the source in all 
segments, although the pulse displays a small variation 
with time in the flrst 300 ks and sharply increases to 
a second almost constant level in the last 200 ks. This 
step-like behaviour is found for most of the spectral 
parameters, with opposite trends. The luminosity of 
the pulse remains remarkably stable with time, in the 
MECS range (1.6-10 keV) despite the variation in the 
relative best-flt parameters. 

— The disk luminosity increases in the pulse and the 
Rin — Tin correlation can be more accurately inter- 
preted using the slim disk model. 

— The luminosity of the corona is lower at the bursts, 
possibly because the corona condenses onto the disk. 

We stress that some of these results are model- 
independent being mainly related to the observed rate. 
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All models that we have investigated show that the lumi- 
nosity below 10 kcV reaches the same level in all runs at 
the peak, whereas at the same time the flux above 20 keV 
decreases respect to the SLT. This behaviour should be 
taken into account in any more detailed modelings of the 
complex emission of this source. 
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Appendix A: Spectral models 

We summarize the results obtained with several combi- 
nations of spectral functions in XSPEC, which did not 
provide acceptable results for all segments. 

i) - Multi-temperature black-body disk -/- power law: 
this is one of the simplest models frequently adopted 
in the spectral ana l ysis of GRS 1915 +105 observations 
( Belloni et allllQQTt Ixaam et al.lll997t ). We obtained ac- 
ceptable reduced only for a few spectra and, in all data 
series, large residuals were found to be apparent above 10 
keV, which cannot be due to instrumental errors because 
they were not eliminated by a change in the intercalibra- 
tion factors within the acceptable range. Adding an expo- 
nential cut-off to the power law improves the fits giving 
narrower reduced distributions, which are marginally 
compatible (^2.5 Gaussian a far) with the expected one 
for the spectra in the SLT-1, SLT-2, and FDT segments 
where an average spectral index of 2.4 and a cut-off at 60 
keV are measured. However, large residuals above 30 keV 
in the P-1 and P-2 spectra, caused by the underestima- 
tion of the flux, produce distributions of the reduced of 
negligible probability of being equivalent to the expected 
one. 

ii) - Multi-temperature black-body disk + Comptonised 
spectrum: the change i n the cut-off power law with a comp- 
tonization spectrum ( Sunvaev fc Titarchuk 1980l ) from 
high temperature electrons (compst) provides the same 
results as the previous model. 

Hi) -Multi-temperature black-body disk + power law 
+ reflection from accre t ion disk: according to previ- 
ous results ( Feroci et al.l Il999l: iMartocchia et al.l l2002l ). 
we added the Compton reflection from an accretion 
disk that is either cold (pex rav) or ionized (pexriv) 
(jMagdziarz fc Zdziarskil 119951 ) to the multi-temperature 
black-body disk plus power-law model. These models 
did not originally include any iron emission line which 
was added with a Gaussian of fixed width (0.1 keV). 
Acceptable fits were obtained for the i^DT interval, where 
the spectral statistics are the lowest, and in the SLT-1 
spectra before the flux increase at 380 ks; the x^ distri- 
bution in these cases is compatible to within 2a with the 
expected one. The two models pexriv and pexrav pro- 
vide compatible best-fit parameters, and in particular the 
ionization parameter in pexriv assumes values compati- 
ble with zero. In the other segments, the model is unable 
to describe the data above 15 keV, where the highest resid- 
uals are observed. 

iv) - Comptonization from bulk-motion matter. 
t he model (BMC in XSPEC) was developed by 



in energy originating from the dynamic bulk motion. The 
co mbination of two BMC co mponents, which was used 
by iTitarchuk fc Seifinal (|2009f ) to fit many RXTE spec- 
tra of GRS 1915+105, is characterized by the black-body 
temperature of the soft photon source, a spectral (energy) 
index, and an illumination parameter related to the frac- 
tion of the bulk-motion flow irradiated by the thermal 
photon source. We left the temperatures of the two BMC 
as free parameters, because fixing them to the value used 
bv ITitarchuk fc Seifina (2009) for the RXTE data (1 keV) 
did not provide acceptable fits to all spectra. An iron line 
with a Gaussian shape was added to the model, together 
with an high energy cut-off for one of the two BMC com- 
ponents. Once again, we found that the best-fit spectral 
distributions, particularly those of the P-1 and P-2 inter- 
vals, were not acceptable, in term of either the x^ values 
or the large residuals. 

The use of an equivalent column density, Nh 
= 5x10^^ cm~^, as the one already considered by 
Titarchuk fc Seifina (l2009l) . provided only a marginal im- 



Titarchuk et al.l (|1997| ) to describe the emission of accret- 
ing matter onto a black hole and takes into account the 
Comptonization of soft photons by a hot gas and the gain 



provement to the fits. The spectra from the SLT-1 and 
FDT intervals could be modeled by a combination of two 
BMCs, the x^ distribution of the spectra were compati- 
ble at the 2a level with the expected ones, and the spec- 
tral parameters of the first BMC are well-defined with an 
approximately stable temperature {kT ~ 1 keV), photon 
index (F = 2.8), and Comptonised fraction (/ = 0.5). The 
second BMC has a flux about 60% fainter than the first 
and a lower temperature (kT ~ 0.5 keV), but the pho- 
ton index and the Comptonised fraction are not confined 
for almost all spectra. Spectra for the SLT-2 intervals are 
marginally compatible with the model, their distribution 
~2.5 Gaussian a being far from that expected, while large 
residuals above 20-30 keV are present in the pulse spectra 
(see Fig. 5). Moreover, the x^ distribution is not compat- 
ible with the expected one. 

v) -Multi-temperature black-body disk + power law + 
black body: To reduce the large deviations in the residu- 
als between 10 and 30 keV obtained with the other com- 
posite phenomenological models, we added to the multi- 
temperature black-body disk plus power law a third spec- 
tral component whose emission is mainly concentrated in 
this limited band. This model cannot be rejected based 
on statistical arguments giving acceptable reduced x^ and 
flat residuals for all spectra in all segments. The results for 
the disk are in close agreement with those relative to the 
"good model" . The power law shows a flux that decreases 
at the peak coherently with the results of the disk+hybrid 
corona model and the spectral index displays small varia- 
tions (<10%) with both time and segment, having a mean 
value of a=2.6. The black-body temperature ranges be- 
tween 3 and 6 keV with the highest values being observed 
during the SLT; its contribution to the source luminos- 
ity is always very low, but becomes significant at ener- 
gies higher than 10 kev, particularly in the pulse, when 
the power law is faint. No correlation with the power- 
law luminosity is observed during the bursts, and we note 
moreover that during these periods it subtends a different 
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spectral range, across which the temperature decreases to 
3 keV. However, a consistent scenario to explain the na- 
ture and t he changes in this compon ent cannot be easily 
envisaged. Titarchuk fc Seifina also found a similar 

component with a temperature of 4-5 keV in some of the 
spectra observed with RXTE and proposed an interpreta- 
tion in terms of a highly gravitationally redshifted anni- 
hilation line. Such a high redshift can be reached only if 
the emission region is at an extremely short distance from 
the Schwarzschild radius, and therefore it seems quite un- 
likely that the flux would exhibit the same modulation of 
the disk. 
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Table 1. Log of the BeppoSAK runs used for the analysis. The columns list the observation codes, the name of the 
runs, the time in seconds from the starting time of tlic first observation (20 October 2000 at UT = 21'* 26"* 55*), the 
exposure time, the rate observed in the energy range used in the analysis, and the number of run bursts. 



Obs. code Run Tstart Exposure Rate N burst 

(s) (s) (ct 











MEGS 


PDS 


HP 




21226001 


A2b 


8994.0 


2191.0 


203.1 


39.6 


80.2 


43 




A3 


14727.0 


2509.5 


203.9 


39.2 


82.6 


53 




A4 


20648.0 


3152.0 


198.0 


36.5 


80.1 


66 




A5 


26511.0 


2809.5 


193.4 


35.1 


79.1 


57 




A6 


31928.5 


3108.0 


191.8 


35.4 


77.2 


64 




A7 


37662.0 


3156.5 


199.2 


38.4 


79.8 


63 




A8b 


43402.0 


2507.0 


197.2 


38.4 


78.7 


62 




A9 


49296.0 


2952.5 


206.1 


41.9 


82.5 


56 


212260011 


B2b 


94999.0 


2197.0 


207.6 


40.6 


_ 


45 




B3 


100740.0 


2638.5 


199.6 


38.0 


— 


55 




B4 


107594.0 


2008.0 


207.6 


40.5 


_ 


38 




B5 


112199.0 


3093.0 


212.2 


42.5 


_ 


61 




B6 


117933.0 


3146.5 


208.3 


40.7 


— 


65 




B7 


123666.5 


3122.0 


199.3 


38.6 


_ 


62 




B8 


129400.0 


3057.5 


203.3 


39.2 


_ 


61 




B9b 


136123.5 


2130.0 


205.8 


41.0 


— 


42 




BIO 


141260.0 


2685.0 


210.1 


41.6 


- 


52 


212260012 


CI 


181003.0 


2198.0 


203.4 


36.5 


83.8 


34 




C2 


186736.5 


2679.0 


206.1 


41.1 


84.0 


49 




C3 


192470.0 


3103.5 


206.6 


39.2 


83.3 


61 




C4 


198203.5 


3121.0 


207.3 


40.8 


83.6 


61 




C5 


203937.0 


3086.5 


210.9 


43.0 


84.6 


59 




C6 


209816.5 


2942.0 


211.6 


40.2 


86.0 


59 




C7 


215404.0 


3067.5 


212.5 


39.7 


88.1 


60 


212260013 


D5b 


267005.0 


2219.0 


215.3 


40.9 


86.3 


43 




D6 


272748.0 


2675.0 


211.3 


39.7 


85.5 


54 




D9 


289938.0 


3117.5 


228.7 


39.7 


96.5 


50 




DIO 


295671.0 


3129.5 


227.8 


39.0 


96.8 


45 




Dll 


301404.5 


2987.5 


232.2 


39.4 


98.0 


39 




D12 


307138.5 


3158.0 


225.8 


40.2 


98.4 


40 




D13 


313230.0 


2730.5 


216.5 


39.1 


91.8 


40 


212260014 


E4 


358737.5 


2666.5 


221.3 


38.9 


91.5 


39 




E5 


364470.5 


3063.0 


212.9 


37.5 


88.7 


47 




E6 


370204.0 


3025.5 


204.7 


39.6 


85.3 


42 




E7 


375937.0 


3112.5 


201.9 


40.8 


83.1 


49 




E8 


381670.5 


3091.0 


212.3 


42.7 


84.8 


59 




ElO 


393136.5 


3062.0 


215.8 


44.2 


84.9 


55 




E12 


405336.0 


2361.5 


235.1 


51.2 


94.7 


36 


212260015 


Fl 


439021.0 


2194.0 


243.0 


52.5 


97.9 


31 




F2 


444735.5 


2644.5 


231.3 


48.9 


92.3 


40 




F3 


450468.0 


3020.0 


231.9 


49.5 


93.6 


45 




F4 


456201.0 


3102.0 


241.8 


51.6 


96.3 


47 




F6 


467752.0 


2979.5 


237.8 


50.6 


95.1 


49 




F7 


473400.5 


3099.0 


234.3 


49.4 


93.0 


51 




F8 


479133.5 


2934.5 


233.6 


49.2 


92.5 


46 




F9 


485200.0 


2744.5 


229.6 


48.7 


91.5 


44 




F17 


530768.0 


2578.0 


244.3 


51.9 


97.8 


36 


212260016 


G3 


547929.0 


3070.5 


243.0 


52.1 




46 




G4 


553723.5 


3026.0 


234.1 


48.4 




50 




G5 


559395.0 


3072.5 


235.6 


48.9 




53 




G6b 


566000.5 


2153.0 


244.1 


53.0 




33 




G7 


571145.0 


2781.5 


230.0 


48.4 




42 
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Fig. 5. Left panels: Data and fitting model plus residuals for spectra in run E7. Right panels: Best-fit models (continuous 
line) with the relative spectral components: disk (dashed line) and Compton reflection (dotted line). 
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Table 2. Average best-fit parameters for the BeppoSAX spectra fit with the "good model" in the three intervals I 
(0-170 ks), II (170-380 ks), and III (400-600 ks). The last section of the table shows the results relative to the two 
RXTE spectra. 



Disk 

(keV) (km) 


Corona 

(keV) 


Nh 


Interval I 


SLT-1 1.17±0.02 54.0±2.1 25.7±1.3 
SLT-2 1.15±0.02 66.2±2.5 36.8±1.5 
P-1 1.51±0.04 50.0±3.7 62.0±3.3 
P-2 2.10±0.07 24.0±1.4 53.6±1.5 
FDT 1.51±0.05 31.3±1.8 24.0±1.6 


25.6±1.3 0.8±0.1 9.0±0.9 
26.4±1.1 0.6±0.1 10.1±1.3 
23.2±1.9 0.4±0.1 9.5±0.7 
18.9±1.1 0.4±0.1 6.8±0.8 
24.9±1.2 0.9±0.2 8.7±0.6 


4.64±0.06 
4.98±0.04 
5.27±0.08 
4.86±0.08 
4.31±0.10 


Iut(-r\-al 11 


SLT-1 1.13±0.05 57.0±4.2 24.7±1.7 
SLT-2 1.11±0.02 70.9±3.1 36.7±1.9 
P-1 1.66±0.13 42.2±7.1 61.5±2.8 
P-2 2.07±0.07 25.3±1.1 58.4±4.8 
FDT 1.46±0.07 33.5±2.3 24.2±2.0 


26.1±1.2 0.8±0.1 8.8±0.8 
26.6±1.4 0.6±0.1 10.1±0.8 
19.6±1.6 0.6±0.1 10.2±0.7 
21.2±1.8 0.3±0.1 7.9±0.5 
25.0±1.0 0.9±0.2 8.3±0.8 


4.6±0.1 
5.0±0.1 
5.2±0.1 
5.0±0.1 
4.4±0.1 


Interval III 


SLT-1 1.21±0.02 57.5±2.1 33.6±0.8 
SLT-2 1.17±0.02 71.2±2.4 45.3±0.9 
P-1 1.35±0.02 65.2±2.0 67.1±1.3 

P-2 1.78±0.03 36.4±1.5 63.3±2.1 
FDT 1.56±0.04 33.3±1.9 31.0±0.8 


22.5±0.5 l.liO.l 12.6±0.5 
23.0±0.5 0.9±0.1 14.7±0.8 
24.3±0.8 0.4±0.1 13.0±0.7 

20.9±1.0 0.4±0.1 8.2±0.8 
22.3±0.7 l.liO.l 10.6±0.6 


4.83±0.07 
5.10±0.05 
5.39±0.05 

5.15±0.05 
4.53±0.09 


RXTE 


SLT 1.34±0.02 60±3 56.2±1.0 
P 1.50±0.03 56±1 74.6±3.1 


20.5±0.3 1.06±0.06 15.9±0.3 
21.2±0.4 0.60±0.09 15.4±0.1 


5.7±0.3 
5.3±0.3 



*Flux in the range 0.01-200 keV in unit of 10"" erg cm"'' s"^ 
**in unit of 10^^ cm"^ 



